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ABSTRACT 

The cross-correlation of Sunyaev-Zel'dovich effect (SZ) and weak-lensing imaging surveys can be used to test 
how well hot baryons trace dark matter in clusters of galaxies. We examine this concept using mock SZ and 
weak-lensing surveys based on the forthcoming AMiBA experiment and generated from a pre-heated cosmo- 
logical Af-body/hydrodynamic simulation. A cross-correlation diagram between matched lensing convergence 
and Compton y maps exhibits butterfly-wing-like structures, corresponding to individual clusters, that encode 
rich information about the distributions of hot gas and dark matter. When the cluster redshift and temperature 
are available the slope of a wing reveals the cluster gas fraction and the width of the wing indicates how badly 
the hot gas traces dark matter. On the basis of simulated data we discuss systematic errors in the projected gas 
fraction estimates that would be obtained from such survey comparisons. 

Subject headings: cosmic microwave background — cosmology: theory — dark matter — galaxies: clusters: 
general — gravitational lensing 



1, INTRODUCTION 

Clusters of galaxies, the largest virialized systems known, 
are key tracers of the matter distribution in the Universe. Clus- 
ters are composed of the hot, diffuse intracluster medium 
(ICM), luminous galaxies (stars), and a dominant mass com- 
ponent in the form of the unseen dark matter. The mass ratio 
of hot gas to dark matter in clusters (the gas fraction) is be- 
lieved to be a fair estimate of the universal baryon fraction 
(White et al. 1993). However, non-gravitational processes 
associated with cluster formation, such as radiative gas cool- 
ing or pre-heating, would break the self-similarities in cluster 
properties and hence cause the gas fraction to acquire some 
mass dependence (Rravtsov et al. 2005; Bialek et al. 2001). 
Weak gravitational lensing on background galaxies probes di- 
rectly the projected cluster mass distribution (Bartelmann & 
Schneider 2001). The thermal Sunyaev-Zel'dovich (SZ) ef- 
fect on the cosmic microwave background spectrum, on the 
other hand, measures the projected thermal electron pres- 
sure (Birkinshaw 1999). Therefore, a combination of weak- 
lensing and SZ measurements will allow us to measure di- 
rectly the cluster gas fractions when the ICM temperature is 
available (Myers et al. 1997; Holder, Carlstrom, & Evrard 
2000). 

Both blind SZ surveys, such as that to be done with AMiBA 
(Ho et al. 2004), and gravitational weak lensing surveys, such 
as the CFHT Legacy Survey 4 , aim to detect massive clusters 
of 10 I4 M Q to study the formation of structure through the 
abundance and properties of the detected clusters. A com- 
parison of the two types of survey should reveal much about 
the different structures of hot gas and dark matter in clusters 
of different mass. Here we investigate this possibility using 
A/-body simulations. 

2. CLUSTER GAS FRACTIONS FROM WEAK-LENSING AND SZ 
OBSERVATIONS 
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The lensing convergence k is essentially the surface mass 
density projected on the sky, 



K — S m /S cr it(Zrf,Zj) ) 
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where £,„ is the surface mass density of a halo acting as grav- 
itational lens of redshift Zd, and £ cl -i t is the lensing critical sur- 
face mass density as a function of lens redshift Zd and source 
redshift z s specified for a given background cosmology (e.g., 
Bartelmann & Schneider 2001). 

The SZ effect is described by the Compton y parameter de- 
fined as a line-of-sight integral of the temperature-weighted 
thermal electron density (Birkinshaw 1999): 
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where n e is the electron number density, T g the gas temper- 
ature, m e the electron mass, crj the Thomson cross section, 
Yj g the surface gas mass density, and X = 0.76 the hydrogen 
abundance. Note that we have assumed isothermality in the 
second equality in Eq. (0. The presence of temperature gra- 
dients in the ICM still remains controversial (see Piffaretti et 
al. 2005 and references therein). In the paper, we assume the 
ICM is isothermal as a zeroth-order approximation; we found 
a temperature variation within the virial radius by a factor two 
at most from our sample of massive halos (see §4). 

In general, the lensing k and the Compton y reconstructed 
from noisy observations are smoothed by a certain window 
function in the map making process. In interferometric obser- 
vations, the synthesized beam (or PSF) is well approximated 
by a Gaussian window. In weak lensing surveys, Gaussian 
smoothing is often adopted to find clusters as local maxima 
in the reconstructed convergence map (e.g., Miyazaki et al. 
2002). We hence express the convergence field smoothed with 
a Gaussian window Wg(0) = exp(-# 2 / 'Oq) / 'tt6 g as kg{0); 0g 
is related to FWHM of the Gaussian window 0f w hm by Og = 
0fwhm/-\/41n(2). Similarly, we define the Gaussian smoothed 
Compton y parameter, yd®)- We then compare the recon- 
structed kg and yc values at a certain point 9. Taking the ratio 
of smoothed kg andyc, we have 
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where f g (9) = £ g .G(<?)/£,„,G(0) is the local gas fraction de- 
fined with the Gaussian smoothed surface mass densities. In 
order to extract information on the cluster gas fraction from 
weak lensing and SZ observations, we thus need additional 
information on the ICM temperature T g , the cluster redshift 
Zd, and the redshift distribution of background galaxies. Such 
information can be in principle available from X-ray observa- 
tions and photometric redshift measurements of background 
galaxies. 

We can define a global estimator for the cluster gas fraction 
based on weak-lensing and SZ observation as 

(ri g ) = ^ iT} * J oc £ crit (z dA ) (f g T s ) (4) 

where r) gl t = rjgiO,) and w,- = w(6j) is the weight for the ith 
pixel Oj. Accordingly, the global gas fraction of a halo can be 
estimated from k and y maps as (f g ) oc (r) g ) T~ l E^ t if one as- 
sumes the isothermality of the ICM. The case with w,- = kg(0,) 
corresponds to a commonly used definition of the gas fraction, 
r/ l) = (y) / (k) oc (E g )/(S m ), as the ratio of the mean surface 
mass densities. However, the mean convergence (k) is ill con- 
strained from weak-lensing observations: Weak lensing k re- 
constructions based solely on image distortions suffer from 
the mass sheet degeneracy (Schneider & Seitz 1995) unless 
additional information such as the magnification bias is taken 
into account (Broadhurst, Taylor & Peacock 1995; Broadhurst 
et al. 2005). On the other hand, interferometric observations 
are insensitive to the DC signal, and hence do not constrain 
the total flux in the field-of-view (FoV). Bolometers also suf- 
fer from a high level of contamination by atmospheric and 
other environmental signals, and hence rely on the differenc- 
ing scheme to extract the sky signals. We therefore consider 
an alternative weight w,- = k g (<?,) that weights high density re- 
gions more strongly than w, = kg(#,). Then we have the gas 
fraction estimator of the form: 

( 2) = J2i K c(0i)yG(0i) = (Ky)_ 

'* ~ £, kg(0,)kg(0,) (k 2 )' ^ ' 

Suppose that k and y satisfy a linear bias relation y{9) = 
an(9). Then errors in the Compton parameter of the form 
y — > y + e (where y 3> |e|) lead to errors in the two es- 
timators of rig of A rfp = c/(k) and Aryj, 2 ' = 6(k)/(k 2 ). 
The difference between these errors lA?^ 1 '! - |A?7^| = 
(l e l/( K )) (Var[K]/(K 2 )) > 0, with equality only if variance 
in k vanishes. Thus the rffi estimator is always better than the 
rfp estimator. Similarly, an error on the convergence of the 
form k — > n + e (where k» |e|) leads to |A7~ (1) | - |A?- (2) | ss 
(a\e\/ (k)) (Var[«:]/(K 2 )) > 0, again showing that the cross- 
correlation based estimator rfp is more robust than ^j, 1 '. 

In practical observations, noise contributions to the estima- 
tor @ must be taken into account. Since the noise proper- 
ties between the k and y maps will not be correlated, the re- 
quired correction is in the denominator of Eq. (|5}, where we 
replace (k 2 ) by (k 2 ) -<t 2 with cr 2 being the noise variance in 
the smoothed convergence field, kg (§3.3). 

3. MOCK OBSERVATIONS 

We use cosmological simulation data by Lin et al. (2004) 
to demonstrate what useful information we can extract from 
the cross correlation between weak-lensing and SZ surveys. 
Specifically, we simulate the forthcoming AMiBA experiment 
as an illustration of such combined surveys. 



3.1. N-body/hydrodynamic simulations 

To make sky maps with realistic SZ and weak lensing sig- 
nals, we use results from preheating cosmological simulations 
of a ACDM model (fl,,, = 0.34, A = 0.66, fl b = 0.044,/z = 
0.66, o"8 = 0.94) in a 100/~ _1 Mpc co-moving box generated 
with an N-body/hydrodynamics code GADGET (Springel et 
al. 2001), which reproduce the observed cluster Mx-Tx and 
Lx-Tx relations at z = (Lin et al. 2004). We construct a 
total of 29 SZ sky maps, each ldeg 2 in size and containing 
1024 2 pixels, by projecting the electron pressure through the 
randomly displaced and oriented simulation boxes, separated 
by 100/T 1 Mpc, along a viewing cone out to the redshift of 
z = 2. Similarly, 29 K-maps are constructed by projecting the 
distance-weighted mass overdensity 5p out to a source plane 
at a fixed redshift, z = z s - 

3.2. AMiBA SZ cluster survey 

As an AMiBA specification, we adopt a close-packed 
hexagonal configuration of 19 x 1 .2m dishes on a single plat- 
form. AMiBA operates at 95 GHz, and this array configura- 
tion yields a synthesized beam of #f w h m ~ 2' (0.6/r'Mpc at 
z = 0.8), which is optimized to detect high-z clusters (Zhang 
et al. 2002). The primary beam has an FoV of FWHM w 1 1 '. 
We assume the following telescope system: bandwidth Av = 
20 GHz, system temperature r sys = 70 K, system efficiency of 
0.7, and dual polarizations. 

AMiBA will operate as a drift scan interferometer to re- 
move ground contamination to high order (Pen et al. 2002 5 ; 
Zhang et al. 2002; Park et al. 2003). We approximate the 
resulting synthesized beam in the AMiBA map by a Gaussian 
with 6V w h m = 2', and we convolve the simulated Compton y 
maps with this Gaussian beam to generate a set of noise-free 
AMiBA yc-maps. The expected rms noise over a sky area of 
fl s with a total integration time of ?j nt is (Pearson et al. 2003), 

=2.0x 10- (^-Y" 2 (J\) "\^) 

V 7 0Ky V 2 0GHz y Vldeg 2 / V240hrs/ 



3.3. Weak lensing cluster survey 

We assume a weak lensing survey with a mean galaxy 
number density of n g = 40 arcmin -2 , an rms amplitude of 
the intrinsic ellipticity distribution of a c = 0.3, and a mean 
galaxy redshift of z s = 1, which are close to typical values 
for a ground based optical imaging survey with a magnitude 
limit of R\{ m 25.5 mag in a sub-arcsecond seeing condition 
(e.g., Hamana, Takada, & Yoshida 2004). The rms noise in a 
Gaussian-smoothed convergence field, kg, is given as 

(e.g., van Waerbeke 2000). We choose the same smoothing 
scale of #f w hm = 2' as for the AMiBA SZ survey and convolve 
the k maps fromiV-body simulations (see §3.1) with the Gaus- 
sian beam to produce a set of noise-free kg maps. Note that 
the adopted value of 6>f w hm is close to the optimal smoothing 
scale of 6* tw h m ~ 1.'7 (or Oq « 1') for an efficient survey of 
massive halos with M ,> 1O 14 /! _1 M found by Hamana et al. 
(2004) based on cosmological A^-body simulations. 

5 http://www.cita.utoronto.ca/~pen/download/Amiba 
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FIG. 1 . — Cross correlation between kq and yc maps as a 2D scatter plot 
histogram obtained from 29 X ldeg 2 noise-free sky maps of weak-lensing and 
SZ surveys for the AMiBA experiment. The kq and vg values are normal- 
ized in units of rms dispersions k* and y„, respectively. Halo peaks with 
and y 6y» are indicated with markers according to their red- 
shifts: z < 0.3 (circle), 0.3 < z < 0.5 (triangle), 0.5 < z < 1 (square). 



4. CLUSTER GAS FRACTIONS FROM CROSS CORRELATION 
BETWEEN WEAK-LENSING AND SZ SIGNALS 

We carried out a statistical analysis on the simulated weak 
lensing kg and SZ Compton yc maps of 29 x ldeg 2 for the 
AMiBA experiment. Figure 1 shows the cross correlation be- 
tween the noise-free kg and yc sky maps of the entire 29deg 2 
as a 2D scatter plot histogram. The scatter plot is presented 
in units of rms dispersions k* and y* of the smoothed noise- 
free fields kg(0) and yc(9), respectively; k, = 1.4x 10~ 2 and 
y* = 2.2 x 10~ 6 . Halo peaks with kg ^ 4.5k* and yc ~ 6y* 
(« 3.3mJy/beam at 95 GHz) are indicated in Fig. 1 according 
to their redshifts. These thresholds effectively select massive 
halos with M 2 oo ~2x l0 H h~ 1 M Q (see Fig. 2 of Umetsu et 
al. 2004); Ma c = (47r/3)A c /9 c r\ is the mass enclosed by r& c 
with the critical density of the universe p c . Figure 1 clearly 
exhibits different capabilities of SZ and weak-lensing in clus- 
ter detections: weak lensing depends strongly on id through 
Edt, while the surface brightness of SZ effect is independent 
of Zd- In particular, gravitational lensing is less efficient for 
high-z lensing halos due to pure geometrical effects. As a 
consequence, high-z halos tend to appear in 1ow-k and high-y 
regions. 

Another important feature in Fig. 1 is a number of butterfiy- 
wing-shaped structures extending from the origin to points 
corresponding to the halo peaks. Each blade of a butterfly 
wing corresponds to an individual cluster. Hereafter, we call 
the n/y cross-correlation diagram the butterfly diagram. For 
illustrative purposes, we show in Fig. 2 a butterfly diagram 
from the noise-free result for the most massive halo, which 
has M 20 o = 7.0 x 10 14 A _1 Af Q at z = 0.142. The SZ signal for 
this halo is more extended than the weak lensing signal in this 
cluster region. Furthermore, we can see a clear disparity be- 
tween the shapes of the contours in the k and y maps. 

The butterfly diagram of a single cluster shows a wing that 
contains rich information on the relative distributions of dark 
matter and hot gas. The tip of a butterfly wing shows the 
surface mass densities of both dark matter and hot gas in the 
cluster core, and the bulk of the wing reveals the degree to 
which the hot gas traces dark matter. The wing width serves 
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FIG. 2. — Same as Fig. 1 but for a 25' x 25' region of a single massive 
halo at z = 0.142 found in the simulation. The halo properties such as M500 
and 7500 are also indicated within the figure. Also plotted are the vq and kq 
maps smoothed with 6>f w i ml = 2' Gaussian in the left and the right subpanels, 
respectively. A circle with radius of #500 is shown in each subpanel. 
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FIG. 3. — Distributions of the ratio fg\< 8500)/ fg(< rsoo) between the 
projected (Eq. [5]) and true gas fractions against the halo mass M200 for three 
redshift subsamples, obtained with the same halo selection as in Fig. 1 but 
from a total survey area of ml 20deg 2 . Each data point is measured from a 
pair of the Gaussian-smoothed noise-free kq and ya maps, while each error- 
bar represents an observational uncertainty (1<t) based on 50 Monte-Carlo 
realizations of random noise in kq and vg- The sample means over the noise- 
free measurements are indicated as dashed lines. 



as an indicator of the difference between projected distribu- 
tions of dark matter and hot gas. Most halos shown in Fig. 1 
have significant wing widths. These widths arise because the 
baryon distribution is closer to axisymmetric than the dark 
matter distribution due to collisional relaxation of the baryons. 
The local slope f] g (0) of a wing in conjunction with the cluster 
redshift and temperature can be an indicator of the local gas 
fraction f g (8) (see Eq. |3))- Note that the estimator rfg' de- 
fined in Eq. (0 can be interpreted as a global slope of a wing 
in the least-square sense. 

We use the halo sample with kg ^ 4.5k» and yc ^ 6y* 
derived from our noise-free sky maps of 29deg 2 to com- 
pare the halo-based cumulative gas fractions f g (< r) = 
M g (r)/M(r) with their projected estimates /] 2) (< 0) oc r] g 2) (< 
9)Tl' 1 'S^ t (zd,Zs) from cross correlation between k and y; 8 = 
r/D^izd) with Da(z) being the angular diameter distance to 
redshift z- We make this comparison inside rsoo, and calculate 
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for each halo a gas fraction using Eqs. @ and (0 within an 
aperture of projected radius #500 = rsoo/D/tfe/)- We assume 
cluster isothermality and take the emission-weighted mean 
temperature Tsoo inside rsoo as the observable X-ray temper- 
ature. The weight w, in Eq. ® is set to zero if k, < or 
y, < to exclude low-signal regions where cross correlation 
measurements are of low significance due to projection effects 
and/or observational noise. Further, we remove those halos 
whose centers are located outside the inner 50' x 50' region 
of a ldeg sky map, in order to avoid systematics in the gas 
fraction estimates due to the boundaries. The effective survey 
area is thus about 20 deg 2 

Figure 3 shows the distribution of the ratio ff\< 
#50o)//g(< 7"50o) obtained from the noise-free maps as a func- 
tion of M200 for three redshift subsamples: (a) z < 0.3 (where 
we find 8 halos); (b) 0.3 < z < 0.5 (16 halos); and (c) 0.5 < 
z < 1 (13 halos). The sample means and their standard errors 
of ff\< 5 oo)//g« r 5 oo) are (a) 0.96±0.03, (b) 1.10±0.04, 
and (c) 0.95 ± 0.05. Note that the quoted errors represent the 
intrinsic scatter of the estimator f ( g 2) from the noise-free sig- 
nal maps. On the whole, the estimator f!P(< O500) for the 
gas fraction is close to unbiased, but there is a slight (~ 2a) 
over-estimation bias of f g in 0.3 < z < 0.5. We found that 
most halos with overestimated values of f g are associated 
with mergers. The shock heating of gas in mergers boosts the 
temperature in the region between the merging halos. How- 
ever, we ascribe the lower temperature T-joo to this gas, based 
on the emission-weighted temperature of the entire halo, and 
so overestimate the value of f g for the system as a whole. 
This overestimate of f g appears to be smaller at higher z 
because such halos have smaller angular sizes and so suf- 
fer larger beam smearing effects. If the 5 interacting/merger 
systems in subsample (b) are rejected, then the sample mean 
fg 2) (< 05oo) /fg{< ^00) = 1.01 ±0.03, showing that the bias is 
eliminated. Such a debias should be possible based on X-ray 
observations of the survey field. 

On the other hand, the projection effect in k due to inter- 
vening matter (Metzler et al. 2001) also affects cluster f g es- 
timates made using weak lensing data. This is more severe 
for higher-z clusters due to their lower geometrical lensing 
efficiency and smaller angular sizes. If an intervening mass 
concentration is projected onto a cluster, then the gas frac- 



tion of that cluster could be under-estimated. Furthermore, 
the gradients in 3D cluster gas-fraction and temperature pro- 
files (Ettori et al. 2004) can also induce a systematic bias in 
projected fj, 2) estimates, especially for spatially-resolved low- 
z clusters. Intrinsic gradients in T g and f g are manifest as a 
curvature of the wing in a butterfly diagram, and they cannot 
be distinguished unless detailed information is available on 
the temperature distribution in the ICM. 

The error bars in Fig. (3) are based on 50 realizations of 
Monte-Carlo simulated noise maps which are linearly added 
to the Gaussian smoothed signal maps of kg and vg- Fol- 
lowing van Waerbeke (2000), a correlated noise map of kq is 
generated for each realization using the noise auto-correlation 
function specified by the smoothing scale #f w hm and the ob- 
servational noise rms a K (Eq. Q). A similar correlated noise 
map for yc uses the noise rms a v (Eq. 0). For a fixed halo 
mass, the random error in /* 2) is smaller for lower-z halos with 
larger angular scales, #500, over which pixel noise is averaged, 
but systematic errors will be dominant for these clusters. 

5. CONCLUSIONS 

By simulating the AMiBA experiment, we have shown that 
the butterfly diagram that emerges from a comparison of the 
SZ and weak-lensing signals contains rich information on 
the distribution of the cluster gas fraction within an individ- 
ual cluster and between different clusters (§4). The cross- 
correlation based estimator fj, 2) (< 0) should be used for a ro- 
bust estimate of the cluster gas fraction. Overall, this estima- 
tor is close to unbiased, but some halos show intrinsic scat- 
ter that exceeds the observational errors (§3) due to various 
sources of systematics (§4). The bias for the the intermediate- 
z subsample of halos can be eliminated by removing interact- 
ing/merger systems. Careful selection of clusters should min- 
imize systematic uncertainties in f g . Detailed information re- 
garding the temperature distribution in clusters, derived from 
X-ray observations, improves the calibration of the gas frac- 
tion estimator, and yields valuable information on the thermo- 
dynamic history of the ICM. 

We acknowledge Kai-Yang Lin and Masahiro Takada for 
useful discussions. 
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